Langmuir waves dispersion in semi-relativistic spinless quantum plasma 
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Based on the Darwin's Hamiltonian the many particle quantum hydrodynamics is considered. Force 
field appearing in the corresponding Euler equation is considered in details, contribution of different terms 
of the Darwin's Hamiltonian in the Euler equation is traced. For example, relativistic correction to particle 
kinetic energy leads to several terms in the Euler equation, these terms have different form, and one of 
them has form looking like a term appearing from the Darwin term. So, two different mechanisms give 
analogous contribution in wave dispersion, which differ by sign of terms. Microscopic analog of the Biot- 
Savart law , called current-current interaction, and describing interaction of moving charges via magnetic 
field, is also included in our description. The semi-relativistic generalization of the quantum Bohm po- 
tential is obtained. Contribution of relativistic effects in the spectrum of plasma collective excitations is 
considered. 



I. INTRODUCTION 

There is fast growing of the interests to the theory of 
the relativistic [Q]]- JH] and semi-relativistic |@] quantum 
plasma. In this paper we develop the many-particle quan- 
tum hydrodynamics (QHD) 10]- Hill in semi-relativistic 
approximation, which allows successive analysis of many- 
particle effects. In this way we are going to discuss be- 
tweenness relation quantum, thermal, semi-relativistic ef- 
fects in the system of many charged particles, as the result 
we will present complete theory including described effects 
for the spinless charged particles. 

Spin caused effects appear in theory in the semi- 
relativistic approximation, however they play significant 
role in non-relativistic physical systems, such as, for ex- 
ample, ferromagnetic. Spin dynamics is also very impor- 
tant in physics of the quantum plasma, where electrons 
and positrons are most widespread objects and their spin 
is the inherent dynamical property. In last decade a lot 
of papers were dedicated to study of spin dynamics in 
quantum plasma, especially by means the quantum hydro- 
dynamics and Vlasov-like kinetic equations. However, it 
is very interesting and important to understand the quan- 
tum many-particle physics appearing from consideration 
of the Darwin's Hamiltonian, which is the spinless ana- 
log of the Breit's Hamiltonian 11211 . The Darwin's Hamil- 
tonian contains both the non-relativistic terms, which de- 
scribe particle kinetic energy and the Coulomb interaction, 
and semi-relativistic terms describing relativistic correc- 
tion to particle kinetic energy (RCKE), which should be 
important at studying of relativistic beams in plasma, in- 
teraction energy of moving charges (which is also called 
current-current interaction and presents itself generalized 
Biot-Savart-Laplace law), and the term proportional to the 
Dirac's delta function called Darwin term or Darwin's in- 
teraction. 

Suggestion were made in Ref. [6] that in some cases 
contribution of the RCKE much smaller than the Darwin 
term. However our studies of semi-relativistic effects in 
the quantum plasma and especially quantities entering the 



Darwin's Hamiltonian, based on the quantum hydrody- 
namics method, shows that the RCKE leads to existence 
of a number of terms in the semi-relativistic Euler equa- 
tion, which have various forms described below. One of 
the terms caused by the RCKE has form close to the only 
term brought by the Darwin term. 

Contribution of the RCKE and current-current inter- 
action in the plasma wave dispersion and the electron 
beam have recently been considered [13] in terms of the 
many-particle quantum hydrodynamics method developed 
in Ref.s 0]- fllll . Another method of derivation of the 
QHD equation for systems of charged spinning particles 
was suggested later in Ref.s 111411 . II 1 511 . Some aspects of 
quantum plasma physics were reviewed in Ref. II 1611 . Be- 
fore we pass to detailed description of the quantum model 
we admit that classic description of plasma theories can 
be found, for example, in Ref.s |fl7ll . lfl8ll . New method 
of derivation of classic hydrodynamics equation from mi- 
croscopic particles motion driven by Newton's equations is 
presented in Ref. lfl9ll . This method does not include using 
of a kinetic equation as intermediate step of the derivation. 

This paper is dedicated to comparison of contribution of 
the RCKE, the Darwin term and the current-current inter- 
action in the Euler equation, obtaining of the explicit form 
of the semi-relativistic pressure tensor and it's influence on 
the dispersion properties of the longitudinal waves. 

Our paper is organized as follows. In Sec. II we dis- 
cuss basic Hamiltonian and compare contribution of differ- 
ent terms. In Sec. Ill In Sec. IV the method of obtaining 
of dispersion equation is described, linearized set of the 
semi-relativistic Euler equation is presented. In Sec. V 
dispersion dependence of quantum semi-relativistic Lang- 
muir waves is calculated and discussed. In Sec. VI brief 
summary of obtained results is presented. 



H. THE MODEL DESCRIPTION 

The equations of quantum hydrodynamic are derived 
from the non- stationary Schrodinger equation for system 
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of N particles: 

ihd t ifj(R,t) = Hip(R,t) 
with Hamiltonian 
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The following designations are used in the equation CO: 
Ci, mi-are the charge and the mass of particle, h-is the 
Planck constant and c is the speed of light, D" = —ihd" — 



the potentials of the external 



electromagnetic field, df = V" is the spatial derivatives, 
Gij = 1 /r,j , - is the Green functions of the Coulomb in- 
teraction, Tij = Ti — r 7 , 
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is the Green function of the current-current interaction, 
ip(R, t)-is psi function of N particle system, R = 
(ri, r N ). Let's consider physical meaning of the terms 
in the Hamiltonian £[]). 

The first term in the right-hand side of the Hamiltonian 
is the kinetic energy, the second is the RCKE, the third 
is the potential energy of the classic charge in the exter- 
nal electric field, the fourth is the quantum contribution in 
the energy of the charge being in the external electric field, 
which is called the Darwin term. All these terms are valid 
for each particle, as they describe kinematic properties and 
interaction with external field. The next terms describe 
inter-particle interactions. First of all, by the fifth term, 
the Coulomb interaction is presented. The next term de- 
scribes quantum contribution in the interaction of charges. 
It is the Darwin's interaction, corresponding to the Darwin 
term. And the last term describes the current-current inter- 
action, which is the microscopic analog of the Biot-Savart 
law. 

The fourth and sixth terms describe the Darwin term. 
The fourth term shows semi-relativistic contribution to the 
force acting from external electric field on charged particle. 
The sixth term presents Darwin term between two particles, 
which can be considered as a semi-relativistic addition to 
the Coulomb interaction. If we have deal with interaction 
of two electrons the Darwin term is 



H D = -7r( — ) S(r t -r ), 
Kmc 1 



(3) 



where r$ and Tj are the coordinates of the two electrons. 
The explicit form of the Darwin term was derived from 
quantum electrodynamics scattering amplitude. Now we 
have to compare the two Darwin terms describing in- 
teraction interaction with external field, the fourth term 
in equation £0, which appearing in the semi-relativistic 
limit of the Dirac's equation 11121 . and inter-particle in- 
teraction, the sixth term in equation £Ql. Admitting what 
Aj(l/|rj — Tj\) = — 47Ti5(r i — r 3 ) and introducing micro- 
scopic electric field caused by, for example, particle j and 
acting on particle i, as E,j = — V^ej/ry). We see that 
sixth term in Hamiltonian ([T]l can be represented as 
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Comparing the fourth term in the Hamiltonian £Q and for- 
mula © we get that these terms coincide if m 2 —> oo. 
However if we consider interaction of two electrons we 
have To; = mj and from (0 we get 



Hr 



eh 2 



4m 2 c ; 



V,E„- 



2„2 1 « ' 



what differs in the two times from the fourth term in (Q}. 
Actually discussed terms should coincide due to the prin- 
ciple of field superposition, so we put additional factor two 
in the fourth term in the Hamiltonian ([TJ, but we will keep 
in mind that we can make another choice and follow to 
the consequence of the Dirac's equation. At discussion of 
wave dispersion we will consider consequences of the both 
choices of the coefficient in the Darwin terms. 

We have included Darwin's interaction in the Hamil- 
tonian, and for short references below we introduce 
new function Gy which is defined as = Gy — 
(fi 2 /4m i c i )S(T i -T j ). 

G^ leads to existence of two terms in the force field in 
the Euler equation. Let's consider how they emerge during 
derivation of the semi-relativistic Euler equation, what we 
make differentiating current j, appearing in the continuity 
equation, with respect to time and using the Schrodinger 
equation. One of them appears due to commutation of Gy 
with the momentum operator^" in the current j. Let's point 
out that the operator pf exists in the current j due to oper- 
ator of the kinetic energy in the Hamiltonian £[]). In the 
self-consistent field approximation this term has following 
form: 



Fg = -e 2 nd a / dr' 
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The self-consistent field approximation allows us introduce 
electric field E caused by the charges of the system, which 
has following explicit form 



-eV / dr' 
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where n is the particle concentration, and satisfy to the 
quasi-electrostatic Maxwell's equations: 



V x E(r,t) 







VE(r,t) = 4vrJ] ea n a (r,t) 



(7) 
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where subindex "a" describe the species of particle: elec- 
trons, ions and positrons. In the case of the one species we 
can rewrite the last formula as 



VE(r,t) = 47ren(r,t). 
Now, force field Fc © take form 
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where the concentration under the space derivative repre- 
sents source of the field, which in the case of the first term 
has come to equation (|9j. So, using equation (|9]) for the 
mentioned concentration we come to 



e 2 h 2 

Fg = enE a —nd a {d 13 E p ). (11) 

Presented here form of the second term corresponds to 
the semi-relativistic contribution in the force acting on the 
charged particle from external electric field obtained from 
the Dirac's equation lfl2Tl . 

The second term caused by Cry appears due to the RCKE 
which is a semi-relativistic term. In the self-consistent field 
approximation it appears as 
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As RCKE has semi-relativistic origin we can write G^ in- 
stead of Gij in this term. The RCKE also gives other terms 
in the force field, all of them are presented below in the 
Euler equation. 

Equations ( fTTT l and (fl~2l are very similar, they have 
two differences. The first difference is distinction in ten- 
sor structure and the second one is the fact that equation 
(fl~2l contains concentration under the spatial derivative, 
whereas in formula (fTTT l concentration contains as an ex- 
ternal multiplier. 



IH. QUANTUM HYDRODYNAMICS EQUATIONS 

To trace contribution of different semi-relativistic terms 
in the wave dispersion we introduce dimensionless coeffi- 
cients 9 and £ equal to 1 and marking different terms. We 
put 9 in the second term of formula dTTb and £ in the for- 
mula dT2b . 

The first equation of the QHD set is the continuity equa- 
tion: 



In that equation a function of current j(r, t) = 
n(r, t)y(r, t) is arisen, where v(r, t)- is the velocity field. 

Second equation of the QHD set is the Euler equation, 
but in the semi-relativistic approximation function j(r, t) 
appeared in the continuity equation in the current of par- 
ticles, however in contrast to non-relativistic case we can 
not call it momentum density, thus the Euler equation is 
the equation of evolution of particles current [20]. This 
equation has form 

mn(d t + v?V p )v a + d P aP = enE a + -e a ^nv f*W 
pfr 2 pfr 2 
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x [mn(r',t)v p {r',t)v^r',t) + P fjl {r',t)}, (14) 

where E and B are the electric and magnetic fields, e is 
the density of thermal energy, including quantum part, an 
analog of the quantum Bohm potential, e^ 1 - is the an- 
tisymmetric symbol (the Levi-Civita symbol), P al3 is the 
tensor of pressure, which is the semi-relativistic general- 
ization of the sum of non-relativistic thermal pressure p al3 
and the quantum Bohm potential T afi , ne is the density of 
the thermal energy. In right-hand side of equation (fT4l a 
force field locates. The force field consists of the Lorenz 
force, and specific quantum semi-relativistic terms, which 
will be discussed below. 7r ai a(r, r', t) will be presented ex- 
plicitly and considered below after analysis of the structure 
of P a/3 . The vector potential of magnetic field appears in 
the form 



A™ t (r,t) = ^J dv>G a p{i 



r')n(r',t)^(r',t), (15) 



d t n + Vj = 0. 
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which gives contribution in the Lorentz force, the sec- 
ond term in right-hand side of equation (fT4l along with 
external magnetic field. Magnetic field caused by cur- 
rents B = V x A mt satisfies to the quasi-magnetostatic 
Maxwell's equation: 

„ 47T 

V x B = — j, (16) 



4 



and 

VB = 0. (17) 
Before discussion of pressure tensor P a/3 we going to 
discuss the physical meaning of terms of the force field 
presented in the right-hand side of the Euler equation (fT4l) . 
It is especially important as some of these terms pre- 
sented for the first time. Thus, the first two terms present 
the density of the Lorentz force. The self-consistent part 
of the Coulomb interaction gives contribution in the first 
term. The second term contains contribution of the current- 
current interaction in the self-consistent field approxima- 
tion. We already should admit that all terms in the Euler 
equation are presented in the self-consistent field approx- 
imation. Actually, only part of the whole contribution of 
the current-current interaction came in the Lorentz force, 
it also leads to several other terms. They are sixth-ninth 
terms of the Euler equation. In fact, the terms seven-nine 
already appear in the classic semi-relativistic hydrodynam- 
ics (detailed discussion of the classic theory will be pre- 
sented in 113D), but in the quantum theory these terms have 
more rich structure. First of all they contain contribution of 
the exchange interactions via quantum correlations, which 
do not considered in this paper, but they naturally appear 
in the many-particle QHD, and we neglect them here, con- 
sidering the self-consistent field approximation only. The 
third term caused by the Darwin's interaction, it includes 
contribution of both the inter-particle Darwin's interaction 
and the Darwin term describing action of external electric 
field. The terms four and five present contribution of the 
RCKE. The fourth term has simple structure, it contains 
divergence V 9 of the tensor which is product of particles 
concentration on 'V a E 13 . In the fifth term, which contains 
a number of terms in square brackets, the first set of them 
is the convolution of E 13 with the tensor which is the cur- 
rent of the particles current j, as a part of this current we 
have the pressure tensor P al3 . As the fifth term of the Eu- 
ler equation has semi-relativistic nature we should consider 
non-relativistic part of P a ^. The second set of terms in the 
fifth term is the product of electric field E a on the energy 
density. The energy density was separated on two parts 
there. First of them is the density of the kinetic energy of a 
local ordered motion. We need to say that ne is the energy 
density which consists of two parts: thermal energy and 
quantum contribution - an analog of the quantum Bohm 
potential. In one particle case we lose contribution of ther- 
mal motion and quantum- thermal terms, and we get quan- 
tum terms arising for non-interacting particles, e gives no 
contribution in considered below problem, therefore we do 
not present it's explicit form. 

Comparing the third and the fourth terms in equation 
(fl4l we find that in the linear approximation at £ = 9 
they reduce each other, or at £ ^ 9 they common contribu- 
tion is less than each taken separately, and sign of common 
contribution depend on relation between £ and 9. Conse- 
quently, simultaneous consideration of the RCKE and Dar- 
win term is necessary. The seventh term appears because 



of the simultaneous account of the Coulomb and current- 
current interactions. It is an analog of the term four, each 
of them appears in the equation due to commutation of cor- 
responding terms in the current j with the Hamiltonian of 
the Coulomb interaction and energy of charge interaction 
with the external electric field. 
Explicit form of the tensor P a ^ is 
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(18) 



where v ia is the velocity of i-th particle, and it is the sum 
of the velocity field v(r, t) and thermal velocity U;, a is the 
amplitude of the wave function ip(R,t) = aexp(iS/H), 
velocity of i-th particle v ia connects with the phase of the 
wave function as v ia = dfS/m. The first term in for- 
mula is the non-relativistic thermal pressure, the second 
term in this formula consist of two parts, the first of them 
is the non-relativistic quantum Bohm potential, the other 
terms present semi-relativistic effects. Neglecting thermal 
velocities in the semi-relativistic terms of the pressure ten- 
sor P a p we get purely quantum semi-relativistic pressure, 
which is the semi-relativistic generalization of the quantum 
Bohm potential T a p, which explicit form for ideal gas is 
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We also drop contribution of the current-current interac- 
tion. In the result we have 

v 2 

Papi*, t) = p a p + T a p jT a p 
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into the system of equations and neglecting by nonlinear 
terms, we obtain a system of linear homogeneous equations 
in partial derivatives with constant coefficients. Passing to 
the following representation for small perturbations Sf 

Sf = f(uj, k)exp(—iujt + ikr) 

yields a homogeneous system of algebraic equations. 

The Euler equation (fl4ll is very complicated, thus we al- 
low ourself to present the algebraic form of linearized Euler 
equation 



d a ^ ■ d^A^/n - ■ d a A\/n\ . 



(20) 



the first two terms have non-relativistic nature. The other 
terms are semi-relativistic, most of them are proportional to 
v 2 (r, t)/c 2 , except of the four last terms. Thermal pressure 
Pa/3 does not depend on interaction, so we can use equa- 



n 



k B TS a ?, 



tion of state for ideal gas, and we write p' 
where kg is the Boltzmann constant, T is the tempera- 
ture, S a/3 - is the Kronecker symbol. When P a/3 stays in 
a semi-relativistic term we should neglect semi-relativistic 
part and consider non-relativistic one (fl~9l l only. 

Here we present explicit form of 7T a p(r, r',t), which is 
the part of the seventh term in the force field in the Euler 
equation 



7r a ^(r,r',i) 
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fid 
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to close the QHD set of equations we should find approxi- 
mate connection between 7r Q ^(r, r', t) and other hydrody- 
namic quantities. Calculating 7r Q ^(r, r',t) for the system 
of independent particles we get ir a p( r , r '> t) = 0- Thus, in 
the first approximation we do not need to account contri- 
bution of 7r Q ^(r, r', t) in the QHD equations. 



-tojmn 5v a + ik a [ 3mv^ e + 
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where G al3 (k) is the Fourier image of the Green function 
of the current-current interaction (O, it's explicit form is 



The last term in the Euler equation (TT4l l give a linear term 
due to the linear part of the quantum Bohm potential, which 
is a part of P^ 1 ', but it is equal to zero because of the tensor 
structure of G Q,3 (k). The last term in equation d2TT) gives 
no contribution in the dispersion of the Langmuir waves. 

The electric field strength is assumed to have a nonzero 
value. Expressing all the quantities entering the system 
of equations in terms of the electric field, we come to the 
equation 



U) 



he 
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IV. DISPERSION EQUATION OF QUANTUM 
SEMI-RELATIVISTIC LANGMUIR WAVES 



h 2 k 2 
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To get semi-relativistic effects in the form of analytic 
simple formulas we consider quantum motion of electrons 
on the background of motionless ions. We consider the 
small perturbation of equilibrium state like 



n, 



n 



+ Sn e 



+ y e 



E = + SE, B = + SB, Sp e = 3mv 2 e 5n , 

where m is the mass of the electron. In equations (fl~3l) . 
(fl4l) and the Maxwell's equations ©, ®, (O and (fTTl) . 
v 2 e is the thermal velocity, for the case of degenerate elec- 
trons v 2 e is the Fermi velocity. Substituting these relations 



where u Le is the Langmuir frequency, oo 2 Le = 47re 2 no/m. 
The first term in the right-hand side of (l22l) consists of three 
parts: the Langmuir frequency, contribution of the Dar- 
win's interaction, proportional to 9, and the third is con- 
tribution of the fourth term of the Euler equation, caused 
by the RCKE. The second group in equation (l22l also con- 
sists of three parts: contribution of the pressure (thermal 
motion or Fermi pressure), the next part is the well-known 
quantum Bohm potential, the last part is the contribution 
of the RCKE via the semi-relativistic part of the tensor of 
pressure. In the next section we will discussed properties 
of the obtained dispersion dependence. 
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V. DISPERSION DEPENDENCE OF LANGMUIR WAVES 

In this section we will discuss contribution of term pro- 
portional to 6 — £ in the dispersion of the Langmuir waves. 
Other term have the same form as presented in the previous 
section. 

Following to the Breit's Hamiltonian we put 6 = 1, and 
according to formula (fT2b we also put £ = 1. Consequently 
we have 6 — £ = 0. In the results we have no contribution 
of the Darwin term in the Langmuir wave dispersion, since 
it compensates by one of the terms caused by RCKE, thus 
we have 

However, we can choose 6 = 1/2 if we follow to the semi- 
relativistic approximation of the Dirac's equation, and we 
have £ = 1, which is only choice for £. In the result we get 

a 2 / h 2 k 2 \ ( , h 2 k 2 h 4 k 4 \ , 2 

(24) 

where we find the term arising in the balance of the RCKE 
and the Darwin term. If we do not consider contribution 
of RCKE, i.e. we put £ = 0, we will get an analog of the 
formula d24l) . but with sign plus instead of minus. 

VII. CONCLUSION 

We gave derivation of the many-particles QHD equations 
for the semi-relativistic system of spinless charged parti- 
cles. We got contribution of the RCKE, the Coulomb, Dar- 
win and current-current interactions in the Euler equation 
is obtained. Comparison of the different terms contribu- 
tions is made. It is showed that simultaneous account of the 
RCKE and Darwin term is necessary, because the RCKE 
gives a number of terms having different structure, and one 
of them has structure of the term describing the Darwin 
term. The RCKE leads to the complex structure of the 
pressure tensor. The semi-relativistic part of the pressure 
tensor contains terms proportional to the rate of the veloc- 
ity field to the square of light, and also includes the several 
terms proportional to fi, 4 /c 2 and contain more higher spa- 
tial derivative than in the non-relativistic quantum Bohm 
pressure. 

Using developed approximation of the QHD equations 
we studied dispersion dependence w(k) of the semi- 
relativistic Langmuir waves. We got contribution of the 
RCKE, which gives two term in w(k), and the Darwin 
term, which gives one term. However, only one term 
emerges in the w(k), since one of the terms caused by the 



RCKE cuts the term caused by the Darwin term, so sum 
of them equal to zero. It is valid if we follow to the Breits 
Hamiltonian in choosing of the coefficient's for the Dar- 
win term. In the other case, this term does not equal to 
zero. We discuss how different choice of the coefficient's 
for the Darwin term reveals in the w(k). 
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